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ABSTRACT: A composite material comprising randomly distributed spherical particles of two 
different isotropic dielectric-magnetic materials is homogenized using the second-order strong- 
property-fluctuation theory in the long-wavelength approximation. Whereas neither of the two 
constituent materials by itself supports planewave propagation with negative phase velocity (NPV), 
the homogenized composite material (HCM) can. The propensity of the HCM to support NPV 
propagation is sensitive to the distributional statistics of the constituent material particles, as 
characterized by a two-point covariance function and its associated correlation length. The scope 
for NPV propagation diminishes as the correlation length increases. 
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1. INTRODUCTION 



The descriptions of electromagnetic planewave propagation traditionally encountered in stan- 
dard textbooks generally involve positive phase velocity (PPV) — that is, the phase velocity casts 
a positive projection onto the time-averaged Poynting vector. On the other hand, there is growing 
recognition of the importance of negative-phase-velocity (NPV) propagation, wherein the phase 
velocity casts a negative projection onto the time-averaged Poynting vector [1, 2]. Of the many 
exotic phenomenons that follow as a consequence of NPV, negative refraction has been the focus 
of particular attention because of its scientific as well as technological significance [3]. 

Manifestations of NPV are not readily observed in naturally occurring homogeneous materials. 
In contrast, artificial metamaterials may be conceptualized — and in some instances physically 
realized — which support NPV propagation. To date, experimental developments with NPV- 
supporting homogeneous metamaterials have been limited to wavelengths larger than in the visible 
regime, with the micromorphology based on elements of complicated shapes [4, 5, 6]. 

In a recent study, we proposed a simple recipe for a NPV-supporting metamaterial, based on the 
homogenization of a random assembly of two different types of spherical particles [7]. The two types 
of constituent particles, type a and type b, are each made of an isotropic, homogeneous dielectric- 
magnetic material, the relative permittivities being denoted by e a ^ and the relative permeabilities 
by fJ- a ,b- Provided that e a ^ and fj, a> t lie within certain parameter ranges, with the real parts of 
e a fe being negative-valued and the real parts of fi a ^ being positive- valued (or vice versa), the bulk 
constituent materials do not support NPV propagation whereas the corresponding homogenized 
composite material (HCM) does. The constitutive parameters of the HCM were estimated using 
the well-established Bruggeman homogenization formalism [8]. The Bruggeman approach has the 
advantages other recent approaches involving NPV-supporting metamaterials [9, 10] that (i) the 
constituent particles are not resonant, and (ii) it is not limited to dilute composites. 

The sizes of the constituent material particles play a significant role in determining whether 
or not HCMs support NPV propagation [9, 10]. Using an extended Bruggeman approach [11], 
we recently reported that increasing the particle size diminishes the scope for NPV propagation 
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in the HCM [12]. The role of the distributional statistics of the constituent material particles in 
NPV-supporting metamaterials has also been highlighted recently [10] — it is this topic that we 
explore herein. 

The Bruggeman homogenization formalism takes into account the distributional statistics of the 
constituent material particles only at the lowest order, via the volume fractions of the constituent 
materials. A different approach is provided by the strong-property-fluctuation theory (SPFT) in 
which a comprehensive description of the distributional statistics of the constituent material par- 
ticles may be accommodated. The provenance of the SPFT lies in wave-propagation studies for 
continuous random mediums [13, 14], but it has more recently gained prominence in the homog- 
enization of particulate composites [15]. Within the SPFT, estimates of the HCM's constitutive 
parameters are calculated as successive iterates to the constitutive parameters of a homogenous 
comparison medium. The iterates are expressed in terms of correlation functions describing the 
spatial distribution of the constituent material particles. In principle, correlation functions of ar- 
bitrarily high order can be incorporated. 

The SPFT is commonly implemented at the second-order level of approximation wherein the 
distribution statistics are stated in terms of a two-point correlation function and its associated cor- 
relation length [16]. The correlation length is taken to be large relative to the constituent material 
particles, but small relative to the electromagnetic wavelengths. The electromagnetic responses of 
constituent material particles within a region of linear dimensions given by the correlation length 
are mutually dependent, whereas the electromagnetic responses of constituent material particles 
separated by distances much greater than the correlation length are assumed to be independent 
[17]. 

In the following sections, the second-order SPFT is applied to investigate the role of corre- 
lation length in predicting whether or not an isotropic dielectric-magnetic HCM supports NPV 
propagation. 

As regards notational matters, Re{z} denotes the real part and Im{z} the imaginary part 
of a complex-valued scalar z. The free-space wavenumber at angular frequency to is written as 



3 



k = LOy/e fi , with e and fi being the permittivity and permeability of free space, respectively. 



2. HOMOGENIZATION 



Let us consider the homogenization of two constituent materials described in the previous 
section. Both constituent materials are composed of electrically small spherical particles, which 
are randomly mixed together. Since materials a and b are assumed to be passive, the principle of 
causality — combined with the implicit time-dependence exp(— iut) — imposes the constraint [18] 

Im {e e } > 



}, (£ = a,b). 

Im {fi e } > ^ 

Let fi £ (0, 1) (£ = a, b) denote the volume fraction of phase £, with f a + fb = 1- 



(1) 



The HCM is an isotropic, homogeneous, dielectric-magnetic material. By application of the 
SPFT, estimates of its relative permittivity e HCM and relative permeability [i HC m are provided by 
the iterative refinement of the corresponding constitutive parameters of a homogeneous comparison 
medium. Furthermore, it transpires that the relative permittivity and relative permeability of the 
comparison medium are identical to those yielded by the Bruggeman homogenization formalism 
[15, 17] — namely, e Br and /i Sr [7]. 

At the second-order level of approximation, the SPFT estimates emerge as [15, 17] 
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+ 2/X Br / ' \fib + 2/x Sr 
Here, k Br = k Qy Je Br [i Br is the wavenumber in the comparison medium, with the sign of the square 
root term selected to ensure that lm{k Br } > 0. The distributional statistics of the constituent 
material particles are taken into account via the term 



T= / t(R)R ex.p(ik Br R) dR, 
Jo 



(4) 



t(R) 



which contains the two-point covariance function t(R). We choose the simple step covariance 
function [19] 

1, R < L 

~ , (5) 
0, R > L 

with correlation length L. Across a range of physically reasonable covariance functions, the SPFT 
estimates of the HCM's constitutive parameters are not particularly sensitive to the form of t(R) 
[20]. 



After utilizing the long-wavelength approximation | k Br L/2ir\ <C 1, (2) and (3) may be combined 
with the covariance function (5) to obtain [17] 
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3. NUMERICAL RESULTS 

Let us now explore how the correlation length L affects the propensity of the chosen type of 
HCM to support NPV propagation. We do so by evaluating (6) and (7) for specific values of e a fi, 
Va,b, fa and L. 

To allow for direct comparison with numerical results based on the Bruggeman homogenization 
formalism reported elsewhere [7, 12], let us select the following values: e a = —6 + 0.9i, fj, a = 
1.5 + 0.2i, €b = —1.5 + i, and fib = 2 + 1.2i. In Figure 1, the relative permittivity e HCM and 
the relative permeability \i HC m of the HCM are plotted as functions of volume fraction f a for 
k L € {0,0.1,0.2}. The constitutive parameters of the HCM calculated for k L = are identical 
to those provided by the Bruggeman homogenization formalism [7]. The graphs are constrained 
such that Re, Im{e HCM } — > Re, Im{e a } and Re, Im{fi HCM } — ► Re,Im{^ a } in the limit f a — ► 1; 
and similarly Re, Im{e HCM } — > Re, Imje;,} and Re, lm{[i HCM } — > Re, Im{^} as f a — > 0. The 
correlation length influences most obviously the imaginary part of e HCM : indeed, lm{e HCM } in 
Figure 1 increases markedly as the correlation length increases. 
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The real and imaginary parts of the normalized wavenumber 

kncM _ , ( Q\ 

r — V ChcmH-hcm {o) 

in the HCM are plotted against volume fraction f a in Figure 2, for k L G {0,0.1,0.2}. The 
sign of the square root term in (8) is chosen to ensure that Im{/c HCM } > 0, as befits a passive 
material. Significantly, the real part of k HCM is negative-valued for mid-range values of f a . It is 
also noteworthy that | lm{k HCM } | 3> | Re{k HCM } | for this homogenization example. 

The parameter 

p _ R e {c^TCm} _|_ RG {pHCAl} /gN 

HCM Im {e HCM } Im {/i flCM } 
is used to determine whether or not the HCM supports NPV propagation [21]; NPV is indicated 

by Phcm < 0. In Figure 3, the NPV parameter p HC M is graphed against volume fraction f a for 

k L G {0,0.1,0.2}. The positive values of p H cu in the limits f a — > and / — > 1 confirm that 

neither of the constituent materials a and b support NPV propagation. In contrast, the HCM 

clearly does support NPV propagation for mid-range values of f a . The range of f a values at which 

the HCM supports NPV propagation decreases as the correlation length increases. 

We explored this issue further in Figure 4, wherein regions of NPV and PPV are mapped in 
relation to Re{e a } G (—6,-1) and Im{e a } G (0,1). The other constituent material parameter 
values are the same as those for Figures 1-3; i.e., p a = 1.5 + 0.2i, = —1.5 + i, and = 2 + l.2i. 
The volume fraction is fixed at f a = 0.3 and k L G {0,0.1,0.2}. At k L = approximately half of 
the mapped e a -space supports NPV propagation, but this proportion decreases as the correlation 
length increases. In particular, NPV propagation is supported only for small values of Im {e a } when 
k L = 0.2. 

4. CONCLUDING REMARKS 

The isotropic dielectric-magnetic HCM arising from a random mixture of two isotropic dielectric- 
magnetic materials — neither of which supports NPV propagation itself — supports NPV prop- 
agation within certain parameter ranges. This conclusion, which had previously been established 
by the Bruggeman [7] and the extended Bruggeman [12] homogenization formalisms, is herein con- 



firmed by the more sophisticated SPFT. In contrast to previous studies involving NPV-supporting 
metamaterials [9, 10], the SPFT-based HCM supports NPV propagation across a wide range of 
volume fraction and the prediction that the HCM supports NPV propagation does not rely upon 
resonant behaviour by the constituent material particles. 

By increasing the correlation length, the scope for NPV is found to diminish, but not disappear. 
In this respect, the effect of increasing the correlation length is similar to the effect of increasing 
the size of the constituent material particles [12]. That the correlation length and the particle size 
give rise to similar effects has been observed elsewhere, in a different context [22]. 
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Figure 1: The real and imaginary parts of the relative permittivity e HCM and the relative per- 
meability Hhcm, of the HCM as estimated using the second-order SPFT in the long-wavelength 
approximation, plotted against volume fraction f a for k L = (solid curve), 0.1 (dashed curve) and 
0.2 (broken dashed curve). Constituent material parameter values: e a = — 6 + 0.9i, \i a = 1.5 + 0.2i, 
€b = —1.5 + i, and fj,b = 2 + 1.2i. 
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Figure 2: As Figure 1 but for the real and imaginary parts of the normalized wavenumber k 
in the HCM. 
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Figure 3: As Figure 1 but for the NPV parameter p HC M hi the HCM. Shading indicates the region 
of NPV (i.e., Phcm < 0). 
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Figure 4: Regions of NPV and PPV, as estimated by the second-order SPFT in the long-wavelength 
approximation, in relation to Re{e a } G (—6,-1) and Im{e a } 6 (0,1). Whereas k L = 0, 0.1, 
and 0.2, the volume fraction is fixed at f a = 0.3. Other constituent material parameter values: 
[A a = 1.5 + 0.2i, €b = —1.5 + i, and fif, = 2 + 1.2i. 
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